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We derive a generalized Ginzburg-Landau (GL) functional near the tricritical point in the (T, H)- 
phase diagram for the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) superconducting state, in 1, 2, and 
3 dimensions. We find that the transition from the normal to the FFLO state is of second order 
in 1 and 2 dimensions, and the order parameter with one-coordinate sine modulation corresponds 
to the lowest energy near the transition line. We also compute the jump of the specific heat and 
describe in the one-dimensional case the transformation of the sine modulation into the soliton- 
lattice state as the magnetic field decreases. In 3 dimensions however, we find that the transition 
into an FFLO state is of first order, and it is impossible to obtain an analytic expression for the 
critical temperature. In this case the generalized GL functional proposed here provides a suitable 
basis for a numerical study of the properties of the FFLO state, and in particular for computing the 
critical temperature, and for describing the transition into a uniform state. 
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I. INTRODUCTION 

As it was shown long time ago by Larkin and Ovchinnikov and by Fulde and Ferrell at low temperatures 
and when the magnetic field is acting on the spin of electrons only, a transition from normal (N) to modulated 
superconducting state (FFLO state) must occur. Due to this non-uniform superconducting state formation the 
paramagnetic limit at T = becomes larger than the usual Chandrasekhar-Clogston limit [0, H p (0) — A ° / - , where 

Ao = 1.76 T c is the superconducting gap at T = 0. The {T,H) phase diagram for 3D superconductors was obtained 
by Saint- James and Sarma (^] assuming that the transition N — ► FFLO is of second order. It happens that the FFLO 
state only appears at T < T* ~ 0.56 T c H and that the temperature-dependence of the critical field H FFLO (F) is 
strongly influenced by the dimensionality of the system. Indeed, for instance in the 3D case H F £ LO {0) — 0.755^-, 

and i?|o ?LO (0) = ^ for 2D superconductors 0, while this field diverges in the one-dimensional case ||. 

Up to now there is no conclusive experimental evidence for the FFLO state formation, except perhaps for what 
concerns UBe^ |(|. However, for this heavy- fermion superconductor the applicability of the standard theory of 
superconductivity is not evident. The main reason for the difficulties in observing such state experimentally resides 
in the fact that the orbital effect is usually more important than the paramagnetic one, and that the actual critical 
field is mainly determined by the orbital effect. However, for heavy-fermion superconductors and low-dimensional 
superconductors (when the field is applied parallel to the planes or chains) the orbital effect can be suppressed, and 
thence we deal with paramagnetically limited critical field. 

The problem of exact structure of the FFLO state is not solved yet even in the framework of the model of pure 
paramagnetic limit, except for the II? case where the superconducting order parameter in the FFLO sate is described 
by the Jacobi elliptic function ||, 0, ||]. For this reason in this paper we concentrate on the description of the 
FFLO phase in the vicinity of the tricritical point (T* , H p (F*)) where the characteristic wavevectors of the FFLO 
state are small compared with the inverse superconducting coherence length . In this case, the appearance of the 
non-uniform state is related with a change of the sign of the coefficient (3 at the gradient term (3(Wtp) 2 in the free 
energy. In the standard Ginzburg-Landau (GL) theory the coefficient (3 is positive, but it happens to be a function of 
the external field acting on the electron spins (paramagnetic effect) and vanishes at the point (T*, H p (F*)), and then 
becomes negative for T < F* . A negative (3 means that the modulated state corresponds to lower energy as compared 
with the uniform one. So, in order to obtain the modulated vector one needs to include a term with second derivative 
in the GL functional. Moreover, in BCS theory the vanishing of the gradient term at (F* , H p (F*)) is simultaneously 
accompanied with that of the coefficient 7 of the fourth-order term H . Due to this particular property, one needs 
to add higher-order terms such as, ip 6 and (yip) 2 ip 2 . Note that some particular cases of generalized GL functional 
were considered in M, [01. 
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In this article we derive a generalized GL functional for 1, 2 and 3-D superconductors which provides a thorough 
description of the FFLO state near the tricritical point. Here, we find that in one and two dimensions the transition 
from the normal to the FFLO superconducting state is of second order while it is of first order in the three-dimensional 
case. In any case the state with simple exponential modulation of the order parameter, i.e. tp ~ \tpo\ e. lq r 1 considered in 
ref. flio| ], does not yield the minimum energy and always corresponds to a second-order phase transition. In the ID and 
2D cases it is most favourable for the FFLO state to appear with a simple sine structure, i.e. tp ~ iposm(q-r). For 2D 
superconductors such structure has lower energy when compared with the "2D lattice" , tp ~ tpo (sin(g • x) + sin(g ■ y)), 
a result that justifies the choice made in ref. of the FFLO structure for layered superconductors with an order 
parameter depending only on one coordinate. For 31? superconductors, however, the situation turns out to be more 
complicated. In this case we find that the transition is of first order for the sine state with one, two, and three- 
dimensional modulation. Our approach yields a good basis for a numerical study of the structure and properties of 
the FFLO state. 

II. GENERALIZED GINZBURG-LANDAU FUNCTIONAL 

In the paramagnetic limit the Hamiltonian of the system can be written in the mean-field approximation, in a 
d-dimcnsional space as follows 



H = d 



, dr J E [*£(r)£Mr) + °H *+(r)4 ff (r)' 
(tp(r) ■ $i"(r)$t 1 (r)+ft.c.) 



(1) 



where a = +1(— 1) if the electron spin is parallel (anti-parallel) to the magnetic field H and H = fisH. The 
superconducting order parameter is ^p(r) — A ($i(r)$_i(r)), where A is the electron-phonon coupling constant. 

Near the transition line the order parameter ip is small, thus using the Gorkov procedure Jl^ | for deriving the 
Ginzburg-Landau functional, we write 

F = J d d r^L -^Y^J d d r K2 G~{-u u , n - r 2 )G + {u v , r 2 - n)^{r x )^ (r 2 ) 

J d d ri, 2 ,3,AG+(u-r 1 )G-(r A -r 2 )Gz(r 3 -r 1 )G+(r 3 -r 2 ) x 

^{ri)ip{r 2 )^ + (r 3 )ip + (r 4 ) H 

where the Green functions here are defined by 

1 1 

G + (ui u ,p) = G {u„,p) 



^2m~ £F 

Being close to the tricritical point of the (T,H)— phase diagram, the spatial modulation of the order parameter is 
small so that we can expand ip( r i) around r\ in Taylor series 

1>(n) ~ i>(rx) + (p*i - r^i) ■ ^) V + - ((T*< - 1*1) • tp + 0(3) (2) 

After long but straightforward calculations, the free energy can be rewritten asQ 

F = a + (3 |cV| 2 + 7 |^| 4 + S \d 2 ^j\ 2 + n |</>| 2 |^| 2 (3) 
+r) [(^+) 2 (9V) 2 + *p 2 {d*P + ) 2 ] + v \ipf 

with the coefficients 



1 We have also added the sixth-order term because the fouth order one turns out to be small near the tricritical point. 
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*r, n s ,ts r^s nN(0)K 3 irN{0)K 5 



in all dimensions, and in ID 



in 2D 



and in 3D 



p= irN(0)V*K 3 _ s = 



1 12 



7r^(0)^x 5 
16 



3 7TiV(0)^ 5 

8 16 



7rJV(0)V^fr 5 
80 



fx = 8r) 



M = 877 



7TiV(0)V^5 



^V(0)Fjif 5 



7TiV(0)V|X 5 



where Vp is the Fermi velocity, N(0) the electron density of states, and 



K n = 2T ■ Re 



= — Rc 

7T™ 



v - 



n > 1 



1 1~L 

uj v = n(2p + 1)T, z = - — i- 



2 2vrT 



or in terms of the psi function ^(x) = ^ lnr(x), K 3 and -K5 can be rewritten as 



2T (-1)' 
(2ttT) 3 2! 

and near the tricritical point, we have 



Re 



(* (2) W) 



#5 = 



2T (-1)" 
(2ttT) 5 4! 



Re(V 4) (.?) 



« = -^)-b-(»(5-«^)-»(5-^) 

= -jV(0) (H '^ o) -Im*'(-- A 



2ttT 



~2 2ttT' 



where Ho is the field corresponding to the second-order transition into the uniform superconducting state, and it is 
given by j| 



T 

In — = Rc 
T 



^2 2ttT j V 



Note also that in the expressions for all coefficients, except a, in the functional (||) we may set H = Ho- 

III. MINIMIZATION OF THE FREE ENERGY 

Now, we proceed to study different solutions for the order parameter in one, two and three dimensions, and see 
which solution minimizes the free energy. We will compute the free energy when the order parameter is an exponential, 
and a sine with one, two and three dimensional modulation. 



A. One-dimensional case 



If we choose the exponential order parameter 



then the free energy reads 
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F oxp = (a + (3q 2 + 6q 4 ) ■ |^o| 2 + (7 + (M - 2r,)g 2 ) • |^ | 4 . 

Analysing the coefficient at \ipo\ 2 we see that the field corresponding to the second-order transition depends on the 
wave vector q, and that the actual field is the maximum one. Accordingly, we get 

9max = ~2^. a = "0 = ^ (4) 

The free energy of such solution then reads 

I? -1 (a-a ) 2 

oxp " 2^iV(0) -*T 3 1 j 

Similarly, for the sine modulation ip(x) = ipQ sin(g.a;), we find the free energy 

sm ttA^O) -tf 3 1 ' 

Note that since -K3 < the denominators in these free energies are positive in the region of existence of the FFLO 
state, which implies that the transition is of second order. On the other hand, we see that the sine modulation 
corresponds to the state with (twice) lower energy than the exponential one in one dimension, and is therefore more 
favourable. 

It is important to emphasize here that this sine solution for the order parameter is just a limit of a more general 
solution on the transition line. Indeed, we show below that as we drift away from the transition line, we need to use 
an exact solution for the order parameter in one dimension . 

Minimizing the free energy (^) leads to following equation for the order parameter (which is assumed to be a real 
function) 

V> (4) + a ip - (3ip" - fx (</>(</>') 2 + V'V) + 2 7 V 3 + 3 v ip 5 = (7) 

where we have redefined 

Now using a similar approach to that of ref . |l3j , we demonstrate that the ( Jacobi) elliptic-sine function 

ip(x,k) = Ak sn(j,k) (8) 

is an adequate solution of eq.(^), A being its amplitude, k the modulus of the elliptic sine and £ is some "effective 
coherence length" that will be determined later on. Note that sn[x, k) — > sin(x), as k — > 0, and this just happens on 
the transition line, and that the amplitude of ip vanishes at the transition. Next, owing to the Jacobi elliptic function 
fl4|| , ip(x,k) satisfies the following equation^] 



Differentiating this equation with respect to x, we obtain 



e ■ r + (k 2 + 1) • v - ^ • ^ = ° ( io ) 

and differentiating the latter equation in turn twice and making use of it, leads to the fourth-order equation, 
,,w , (fc 2 + i) ,„ 12 ,,, ,,,2, ,2 ,,a , 12 ,5 6(fc 2 + l) 3 



2 The prime here stands for the derivative with respect to the space coordinate x. 
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Multiplying (Q) by ip and (10) by ip , and adding the results leads to 



1 (^r + *vo - tAi • ^ 5 + -^-^^=0 (12) 



£2 A 2 / r r / ^4 A 4 T ^4 A 2 v ^ 4 

and we also rewrite (KOI) as follows 



^ + <i=+I>!.,_ »!*£>., _o (13) 

Now, the point is to obtain an equation for ^ that could be identified with eq.(0). For this purpose, we add 
eq.fjll]) to the results of multiplying respectively equations (|l3|) and (|l^) by arbitrary coefficients A and B. Then by 
identifying the coefficients in the resulting equation with those in (Q), we obtain the following system of 5 equations 
for the 5 parameters A, £, k, A, and B: 

^_ A{e + lf~Bk^ ~_ (1 + A)(P + 1) (w) 



7 : 



' ' e 

(B — A — 3)(fc 2 + 1) 



£4 A 2 

12 -B „ 4-B 



On the other hand, in one dimension the coefficients a, (3, 7, /I and v are obtained from eq. ([}]), 



7 : 



1 6(^1 - X?) ~ _ 4K 3 

4^3 



10 „ 



In addition, equating the ratio ft- in (|lj) with that in ( |l5| ) we infer that B = 2 and Vp = £ 2 A 2 . 

Moreover, as explained above, on the line of second-order transition the solution must be sin(x), i.e. k — > 0, 
whereupon the coefficients in ( |l4[ ) become 

„ A 1+A „ B-A-3 , . 

« = ^4« = ^2—' 7 ^ ^ A 2 > ( 16 ) 

12-5 ^_ 4-B 

/X_ ^2 A 2 ' V ~ "^4^4" 

Then, using the condition for the second-order transition, 

2 

we find that A = 1. 

At this point we wish to stress that in two and three dimensions one can readily check that the elliptic sine is no 
longer a good solution that transforms into sine on the transition line. Indeed, on this line A — 1, so that the identity 

7 = infered from eq.([l6|) is compatible with eq.(|l5|) that is valid in one dimension. However, this is not true in 

other dimensions since in this identity the right hand side depends on the dimension of space whereas the left hand 
side doesn't, see (||). 

Away from the transition line defined by eq.(p^), the parameter A obtained from ( |l4| ) can be written in terms of 
the "external field" parameter h — -S^ as follows 



5 



Of course, A(k, h) — > 1 as k — > 0, since then h —> j due to eg. (p"7|) . 

Note that in this case we have in fact a family of solutions parametrized by k. Indeed, using eqs.{ fl4| ), jl^ ) we can 
express the parameters A and £ in terms of the modulus k, bearing in mind that B = 2 and Vp = £ 2 A 2 , see above. 

Therefore, in one dimension we have found a solution to the system of equations (141) for all parameters entering in 
the order parameter ip(x,k) described by the general Ginzburg-Landau functional (3). Now, let us compute the free 
energy of this non-uniform order parameter and compare it with that of the uniform state. We shall also study the 
transition between these two states. For this purpose, we insert the solution (j^) into the Ginzburg-Landau functional 
(||) and integrate over a period of the Jacobi elliptic sine, namely AK{k), K(k) being the complete elliptic integral of 
first kind [ajj , thus we obtain the following expression 



F NU (k,h) 



C 



where we have defined 



-k 4 



P 2 + i)(A+i)r 



[10 h - (fc 2 + 1)(11 + A) h + 14fc 2 h] 



(19) 



In{k) = 



dx sin™ x 



K ( k ) J \fl - k 2 sin 2 x 

n v 



or in terms of K(k) and the complete elliptic integral of second kind E(k) , 

k 2 2(1 + k 2 ) E{k) 
K(k) 



J 1 (\ E{k) \ J 2 

h -¥y--K(k))' h -~ 



h = 



3fc 4 

4(/, 2 + l) 7 3 / E(k)\ 



5k 2 



5fc 4 



K(k) J 



The coefficient C = S Vp (—/3) 3 is positive since (3 < 0. 

Note that the most advantageous feature of the free energy (19) is that it depends only on the parameter k once the 
external field parameter h is fixed. Then, upon minimizing this free energy we determine the modulus k and thereby 
we obtain the parameters £ and A from eq.(14), 



C 2 = (l + A(k,h))(k 2 + l) 



A = 



/3 e 

On the other hand, the free energy for the uniform order parameter reduces to 

F v = 5- (atp 2 +7?/> 4 + V^) 

where the coefficient 7 is negative, and hence the transition is of first order here. 
Using the same notation as above, the free energy Fjj becomes 



Fu(h) 
C 



-1 



.1 



(l + Vl-6/i)P-2/i)-/i 



(20) 



which depends on the sole parameter h. 

It is seen in figure 1 that the free energies (|l^) and ( p0| ) converge when the external field h becomes equal to the 
critical value he ~ 0.0925925. The latter marks the transition from the non-uniform to uniform state. Next, we show 
in figure 2 the change of the form and period of the order parameter from sine to soliton lattice as we lower the 
external field. When h < h c , the period of the order parameter is infinite and the superconducting state becomes 
uniform, see figure 3. 

Before considering the two and three dimensional cases, let us first compute the jump in the specific heat at the 
second-order transition in one dimension. In this case the free energy for the sine solution is given in eq.(^), and the 
corresponding jump in the specific heat reads 



AC = 4N(0)T C x 



(Imfr'(i 



1 2-kT> 



dTi.0 
dT 



Re *(2)(i_^ 



(21) 



The characteristic feature of such behaviour is the divergence of the jump at the tricritical point. 
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B. Two-dimensional case 



1. Exponential 

As above, we start by considering the exponential order parameter 

V(r) = 4> e iq - r 

for which the free energy reads 

F = (a + Pq 2 + 8q 4 ) ■ |^ | 2 + ( 7 + (M - 2 V )q 2 ) ■ |V>o| 4 
and upon minimizing over q 2 and I'i/'ol 2 ; this becomes 

1 (a - a ) 2 



For one-component sine 



we obtain 



and finally for the two-component sine 



we get 



Therefore, since K3 < 0, we infer that 



Foxp ~ nN(0) K 3 • (22) 
ip(r) = ipo sin(g.x) 



F - 6 ( a - a °) (9 o) 
" m - ~ nN(0) K 3 { 6} 



ip(r) = ipo {sm(q.x) + sin(g.y)) 



4 (a-a ) 2 
ttN(Q) K 3 



(24) 



This means that the phase with one-component sine modulation is the most stable in two dimensions. Here again 
we find that the coefficients of the quartic term in the free energy are positive, which implies that the transition into 
the FFLO state is of second-order. Our result can justify the choice made in ref. [jn] of the structure depending only 
on one coordinate for the numerical analysis of the FFLO phase in two dimensions. 

Finally, as in one dimension we obtain the jump in the specific heat 

AC = 24iV(0)T c x 1 {2 - \** T > dT ' (25) 
Re *( 2 )(!-^) 

which diverges at the tricritical point. 

C. Three-dimensional case 

The free energy of the exponential order parameter tp(r) — ipo e iq ' r , reads 

F = (a + [3q 2 + Sq 4 ) ■ IVol 2 + (7 + {a - 2v)q 2 ) ■ l^Pof ■ (26) 

Then using the coefficients in eq. we see that the coefficient of the quartic term in ( p^ ) is positive at q 2 = q$ = — , 
thus indicating that the transition is of second order. 
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On the contrary, for the 1,2, and 3 component sine we find that this coefficient is negative and consequently the 
transition is of first order. But, in order to compute the critical temperature it is necessary to add the sixth-order term 
vip 6 . The transition occurs into a structure with well-developped harmonics and to find the corresponding critical 
temperature and the structure of the non-uniform state numerical calculations are needed. 

The functional (Q) provides a good basis for such calculations. Apparently the 3D-lattice state is more stable. 
Really, one could also study the change in the order of the transition as we drift toward the T — point. In fact, it 
was argued in |1| that as T — ► the transition into the sine state becomes of second order, whilst the results of ref. 
jl| showed that at this point the transition into the 3D-lattice state remains of first order, and thus corresponds to a 
higher critical field and lower energy as compared with the sine state with one-dimensional modulation. 

We may speculate that the 3-D-lattice is also more favourable near the tricritical point and along the whole transition 
line. 



IV. CONCLUSION 



We have constructed a generalized Ginzburg-Landau functional for the FFLO superconductors near the tricritical 
point in the (T, ff)-phase diagram, in 1, 2, and 3 dimensions. Using this functional we have shown that the state with 
exponential modulation of the order parameter (jp ~ e lqx ) is always unfavourable and that the transition from the 
normal to the FFLO state is of second order in 1 and 2 dimensions, where the order parameter has a one-component 
sine modulation near the transition line. In one dimension we have shown that upon lowering the external field, the 
order-parameter modulation changes from sin(a;) on the transition line into elliptic sine sn(x), and finally transfers 
to a uniform state. The proposed generalized Ginzburg-Landau functional for ID superconductors can be directly 
applied to the description of the spin-Peierls systems in magnetic field near the corresponding tricritical point [|l6) . 

In 3 dimensions we have shown that the transition is of first order, and contrary to 1 and 2 dimensions, one has to 
take into account higher harmonics of the order parameter to calculate the critical temperature. 
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Figure captions 



• Figure 1: Plot of free energy as a function of the "external field" h = = S Vp (— /3) 3 . The solid line 

represents the free energy of the uniform state. The line in crosses represents the free energy of the non-uniform 
soliton-lattice state. These two curves meet at h c = 0.0925925. 



Figure 2 : Spatial dependence of the superconducting order parameter in the non-uniform ID phase for different 
fields h = ~k The dotted line indicates the solution at h = 0.165, and the solid line represents the solution 
when h = 0.0927, i.e. close to the critical value h c ~ 0.0925925. The dimensionless space coordinate x is defined 



byl = I ' (V^Vf) , and k) = V(1+ ^ } )(fc2+1) - 



• Figure 3 : Variation of the period L of the non- uniform phase as a function of the external field h. The period 
diverges as the external field approaches the critical value h c . The dimensionless period L plotted here is defined 
as L = L ■ ( y/—jVp J , and for the soliton-lattice phase this is L = AK(k) ■ + A(k, h))(k 2 + 1), see text. 
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